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I. INTRODUCTION 

Direct photon production in p + p collisions has long 
been regarded as a fundamental observable [Ij. At lead- 
ing order (LO) in perturbative QCD photons are pro- 
duced directly from the hard scattering of partons and 
are hence independent of non-perturbative effects from 
hadronization. The LO diagrams, shown in Fig. [TJ 
arise from two parton scattering processes: quark-gluon 
Compton scattering and quark anti-quark annihilation. 
Figured!^ upper) shows the fractional contribution of lead- 
ing order parton scattering processes for ^/s = 200 GeV 
p + p collisions at midrapidity (|r/| < 0.35) using the 
CTEQ6 parton distribution functions [2l| . Compton scat- 
tering dominates due to the abundance of giuons relative 
to anti-quarks. Due to the preponderance of giuons in the 
initial state, direct photons have historically been used 
to probe the gluon distribution of the proton d, Q ■ For 
comparison. Fig. (^flower) shows the contributing pro- 
cesses to leading order dijet production using a hadron 
trigger, in this case a 7r°. Their relative contributions 
have a non-trivial pt dependence, in large part a con- 
sequence of triggering on a jet fragment^. In the final 
state, to the extent that leading order Compton scat- 
tering dominates, the direct photon is likely to oppose 
a quark jet. Moreover, the transverse momentum of the 
recoil parton is exactly nominally balanced by that of the 
photon, a feature often exploited to determine the energy 
scale in jet reconstruction [6| . 

Beyond LO photons may be produced by 
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bremsstrahlung from a quark. Although the rate 
for hard photon radiation is calculable at NLO, photons 
must be considered part of the jet below an arbitrary 
fragmentation scale. The contribution of these frag- 
mentation photons depends on the non-perturbative 
parton-to-photon fragmentation functions which are 
poorly constrained relative to the fragmentation func- 
tions into hadrons 0, NLO effects spoil the exact 
transverse momentum balance between the photon and 
the recoil jet which holds at LO. Often the contribution 
from fragmentation photons is suppressed by applying 
an isolation criterion to the photon sample |9l4l3l|. 
Typically, the total energy from hadron production in a 
cone around the photon is required to be small compared 
to that of the photon. The cross section of isolated 
photons may then be compared to NLO calculations 
with the same isolation criterion applied. 

At LO, a pair of hard-scattered partons emerges ex- 
actly back-to-back. However, due to the finite size of the 
proton, each of the colliding partons has a small trans- 
verse momentum on the order of 300 MeV |l3|- At NLO, 
an additional transverse momentum component arises 
from emission of a parton in the initial state. Effects that 
give rise to a non-zero p?^" , defined as the vector sum of 
the outgoing parton transverse momenta, are collectively 
referred to as the kr effect, kx being defined as the trans- 
verse momentum per parton (Ifcrl = IPt''^"^I/>/2). The 
component of fcy in the direction transverse to the out- 
going parton pair causes them to be acoplanar while the 
component along the axis of the outgoing parton pair 
imparts to them a momentum imbalance. 

For observables such as Drell-Yan production, NLO 
calculations are insufficient to describe the magnitude of 
the kx effect, since multiple soft gluon emission gives 
an additional contribution requiring a resummation of 
the perturbative series [ll]. Measurements of (p^*^") 
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FIG. 1: Leading order Feynman diagrams contributing to direct photon production, 
quark-gluon Compton scattering and t and u channel quark anti-quark annihilation. 
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FIG. 2: The fractional contribution of parton scattering pro- 
cesses to (upper) inclusive direct photon and (lower) Tr° pro- 
duction at LO in p + p collisions at -^5 = 200 GeV. The 
processes q + q ^ g + g and g + g ^ 1 + q contribute to tt" 
triggered events at the level of approximately 1% and 0.01%, 
respectively, and are not shown. 



show comparable values for Drell-Yan, dijet and diphoton 
events many of which are compiled in Fig. 30 of [l6l |. 

It was argued in [13] that NLO pQCD is insufficient to 
describe the world data on direct photon cross sections. 
This claim was contested in [T^ in which it was argued 
that that inconsistencies amongst the data sets may in- 



stead be responsible for discrepancies between data and 
NLO. Although progress has been made, most notably in 
the form of the joint- resummation approach p^ . there is 
not yet a resummed theory which successfully describes 
all the data. As a consequence, direct photon data are 
not used in most determinations of the parton distribu- 
tion functions. 

Direct photons are also an important observable in 
heavy- ion collisions [23j. A wealth of evidence indi- 
cates that a hot, dense state of matter characterized by 
partonic degrees of freedom is produced in these colli- 
sions [m . Hard-scattered partons are believed to rapidly 
lose their energy to this QCD medium based on the ob- 
servation of the suppression of high-pT hadrons [H, |23[ , 
a phenomenon known as jet quenching. At LO, direct 
photons play the role of the hard-scattered parton (see 
Fig.[T]) with the distinction that, as color-neutral objects, 
they do not interact strongly with the hot QCD matter. 
The absence of any nuclear modification for direct pho- 
tons \2A] therefore acts as a control measurement for the 
jet quenching phenomenon and constrains possible con- 
tributions from novel sources of induced photon produc- 
tion predicted to arise from the interaction of partons 
with the QCD medium 25, 26]. 

Just as photons may be used to determine the jet en- 
ergy scale in p + p collisions, they may be used to es- 
timate, on an event- by- event basis, the initial energy of 
the opposite-side parton in heavy-ion collisions, an idea 
first proposed in [23|. The distribution of hadrons in 
the away-side jet reflects the so-called medium modified 
fragmentation function which is the product of the jet 
fragmentation in the dense QCD environment. Devia- 
tions from vacuum jet fragmentation, as observed in p+p 
collisions, hence should enable tomographic studies of 
the medium using the energy loss of the away-side to 
probe the density profile of the medium. A first, albeit 
statistics-limited measurement, of direct photon-hadron 
correlations was presented in [28] . For a recent review of 
medium-modified fragmentation functions the reader is 
referred to ^2^. 

The interpretation of direct photon triggered correla- 
tions in heavy-ion collisions necessitates detailed mea- 
surements of such correlations in p+p collisions. The mo- 
mentum balance between the photon and the opposite- 
side parton is spoiled due the fcr effect and by the con- 
tribution of fragmentation photons in the direct photon 



5 



sample. The present work endeavors to study such ef- 
fects. The remainder of the article is organized as fol- 
lows. In section II we describe elements of the appara- 
tus relevant to the measurement of photon-hadron cor- 
relations. In section III the methodology of extracting 
direct photon correlations from the background of de- 
cay photon-hadron correlations is detailed. Section IV 
presents results on tt'^ and direct photon triggered cor- 
relations. Finally, section V interprets the results at the 
partonic level using a simple of LO pQCD calculation 
coupled with phenomenologically-motivated kr smear- 
ing. 



II. DETECTOR DESCRIPTION AND 
PARTICLE IDENTIFICATION 

A. The PHENIX Detector 

The PHENIX detector, described in [lOl, is weh suited 
for jet correlation correlations between photons and 
hadrons. The central arms of the detector are nearly 
back-to-back in azimuth (offset by 22.5°), each subtend- 
ing 90° and covering 0.7 units of pseudorapidity around 
midrapidity. Each arm contains charged particle tracking 
chambers and electromagnetic calorimeters (sij . 

The electromagnetic calorimeters (EMCal) [s^l con- 
sist of two types of detectors, six sectors of lead- 
scintillator (PbSc) sampling calorimeters and two of lead- 
glass (PbGl) Cerenkov calorimeters, which measure elec- 
tromagnetic showers with intrinsic resolution ge/E = 
2.1%®8.1%//Eand0.8%®5.9%//E, respectively. The 
fine segmentation of the EMCal (Ary x Ac/) ~ 0.01 x 0.01 
for PbSc and - 0.008 x 0.008 for PbGl) allows for ttO 
reconstruction in the diphoton decay channel out to pt 
> 20 GeV/c. In order to measure direct photons over 
the range 5 < pt < 15 GeV/c, 7r° and must be recon- 
structed over a larger range of 4 < pT < 17 GeV/c to 
account for decay feed-down and detector resolution as 
described below. Direct photon, tt*^ and 77 cross section 
measurements in p + p collisions are described in [ssl - lssj . 
Photon candidates of very high purity (> 98% for ener- 
gies > 5 GeV) are selected from EMCal clusters with the 
use of cluster shower shape and charge particle veto cuts. 
At large photon px (» 10 GeV/c in the PbSc), clusters 
from tt'^ photon pairs start to overlap. Nearly all of such 
merged clusters, as well as other sources of hadron con- 
tamination, have an anomalous shower shape, and thus 
are removed from the analysis. 

Charged hadrons are detected with the PHENIX track- 
ing system |36| which employs a drift chamber in each 
arm, spanning a radial distance of 2.0-2.4 m from the 
beam axis, and a set of pixel pad chambers (PCI) sit- 
uated directly behind them. The momentum resolution 
was determined to be 5p/p = 0.7% ® 1.0%p where p is 
measured in GeV/c. Secondary tracks from decays and 
conversions are suppressed by matching tracks to hits in a 
second pad chamber (PCS) at a radial distance of ~ 5.0 



m. Track projections to the EMCal plane are used to 
veto photon candidates resulting from charged hadrons 
that shower in the EMCal. 

The data used in this analysis consists of approxi- 
mately 533 million photon-triggered events from the 2005 
and 2006 p+p data sets. The total recorded integrated 
luminosities during these runs were 3.8 (2005) and 10.7 
(2006) pb^^, respectively. Events were obtained with 
an EMCal-based photon trigger, described in jsj, which 
was over 90% efficient for events with a photon or n° in 
the range of energy used in this analysis. 

B. tt'' and rj Reconstruction 

The background for the present analysis consists of cor- 
related decay photon-hadron pairs. In order to measure 
this background, tt'^ and 77 mesons are reconstructed in 
the 2-7 channel, which together are responsible for ap- 
proximately 95% of the decay photons. The invariant 
mass windows for tt" and 77 mesons are 120-160 and 530- 
580 MeV/c^, respectively, as shown in Fig. [3] The rate 
of combinatorial photon pairs is reduced by only consid- 
ering photons of energy > 1 GeV. 




Invariant Mass [MeV/c^] 

FIG. 3: Di-photon invariant mass distribution for pairs of 
5 < Pt < 15 GeV/c demonstrating tt" and rj reconstruction. 
The signal and side-band regions are indicated. Offline event 
filtering cuts are responsible for the features at 0.105, 0.185, 
0.4 and 0.7 GeV/c^ 

As discussed in the next section, jet correlations of 
decay photon trigger particles with charged hadron part- 
ner particles are estimated from those of their parent 
mesons. The quantity of interest is the per-trigger yield, 
Y . For the tt", for example, the per-trigger yield, Yt^, is 
the number of Tr'^-hadron pairs divided by the number 
of 7r° triggers. For the given photon selection, the effect 
of combinatorial pairs on the measured per-trigger yield 
was evaluated with PYTHIA to be < 2%. The smallness 
of the effect is due both to the small combinatorial rate 
itself as well as the similarity of correlations from falsely 
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matched photons, most of which are themselves from tt" 
decay, to the true n° correlations. 

In contrast, reconstruction of the rj meson has a much 
smaller signal-to-background of 1.4-1.6, depending on 
the pt selection. In this case, the per-trigger yield of 
the combinatorial photon pairs is estimated from photon 
pairs with invariant mass in the side-band ranges of 400- 
460 and 640-700 MeV/c^, beyond 3cr of the rj peak. The 
procedure for subtracting this combinatorial contribution 
is discussed in [l^ and gives rise to a 10% systematic un- 
certainty on the T] per-trigger yields. 



III. SUBTRACTION OF DECAY 7-HADRON 
CORRELATIONS 

A. Statistical Subtraction 

A direct photon is defined here to be any photon not 
from a decay process including those from parton to pho- 
ton fragmentation. The relative contribution of direct 
and decay photons is expressed in terms of the quantity 
, where R-y = (number of inclusive photons) / (number 
of decay photons). This quantity has been determined in 
the course of the measurement of the direct photon cross 
section in y/s = 200 GeV p+p collisions presented in [33l | 
and its values are given in Table |T]^. In J2§\ the yield of 
charged hadrons per direct photon trigger was estimated 
by a statistical subtraction according to 



^^ircct — j-> (-^7 -inclusive -^ccay)! (1) 

lij — i 

where the trigger particles in the per-trigger yields are 
inclusive photons, decay photons and direct photons as 
indicated. 

The contribution of decay photons to the inclusive 
photon-hadron correlations is determined by weighting 
the measured 7r° (F^r) and rj (F^) correlations by their 
probability to contribute to a given Idecay bin. For a 
given 7r° distribution the number of decay photons from 
7r° in the decay photon bin of (a < < b) is given by 

^7(7r)(a <p^ <b) = 

/ e^{pl,PT) ■ VApI,Pt) ■ ^^\pt) ■ NAPt) dp^T- (2) 

J a 

Vtt{PttPt) is decay probability density for a tt" of p^^ 
to decay into a photon of . e-y and e^r are the single de- 
cay photon and tt'' reconstruction efficiency, respectively. 
The dependence of on p^ is due to cluster merging at 



^ Note that the values of R-y are not corrected for efficiency losses 
due to cluster merging. 



high Pt- For brevity, we define the efficiency corrected 
decay probability density: 



PApI,Pt) ^ ^,{pI,Pt) ■ Wt,Pt) ■ e~\p^). (3) 

For a finite 7r° sample the integral in Equation [5] becomes 
a sum and the per-trigger yield of decay photons from 
is calculated according to 

^ j:N^pipi,PT)N.-h 

'^''^ T.n^P{pI,P^t)N. ' ^ ' 

where the pt limits of the bins have been made implicit 
for brevity. For a perfect detector, P is calculable ana- 
lytically. A Monte-Carlo (MC) generator implements the 
PHENIX acceptance and uses Gaussian smearing func- 
tions to simulate detector resolution according to the 
known EMCal energy and position resolution. Occu- 
pancy effects give rise to an additional smearing of the 
TT*^ and rj invariant masses. This effect is included in the 
MC by tuning the resolution parameters to match the 7r° 
peak widths observed in data. The uncertainty on the 
decay photon mapping procedure, including the effect of 
combinatorial photon pairs in the t:^ matching window, 
was evaluated in PYTHIA to be 2%. This procedure is 
described in more detail in [28j. 

Once the decay photon correlations for tt*' and rj trig- 
gers have been obtained, they are combined according 
to 



i^dccay = (l/fJ^M) Y^(^) + (l - I/'^tw) Y^i-n)- (5) 

The quantity (5^(^) is the ratio of the total number of de- 
cay photons to the number of decay photons from 7r° . Its 
value was estimated to be 1.24 ± 0.05 33]. Inefficiencies 
in the detection of photons form i:^ decay in increase 
the ^^(Tr) slightly, giving a value of 1.28 for decay pho- 
tons in the range 12-15 GeV/c. Equation [S] effectively 
assigns the same per trigger yield to the heavier mesons 
{u!,rji,4>,...) as for rj triggers. This assumption was stud- 
ied in PYTHIA and found to influence Idccay at the level 
of 2%. The total systematic uncertainty on the decay 
photon associated yields contains contributions from the 
?7 sideband subtraction, the value of (5-y(7r), the effect of 
hadrons heavier than rj and the MC decay photon map- 
ping procedure, in approximately equal parts. 

The per-trigger yields are corrected for the associated 
charged hadron efficiency using a GEANT simulation of 
PHENIX detector. The quoted yields correspond to a 
detector with full azimuthal acceptance and \rj\ < 0.35 
coverage. No correction is applied for the A77 acceptance 
of pairs. A pr-independent uncertainty of 8% was as- 
signed to the charged hadron efficiency. 
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B. Decay Photon Tagging 

The statistical and systematic uncertainties of the di- 
rect photon correlation measurement may be improved 
by event-by-event tagging of decay photons. It is not 
possible, however, to remove all decay photons due the 
finite acceptance and efficiency for parent meson recon- 
struction. The remainder must be subtracted according 
to the statistical method described in the preceding sec- 
tion. Photons are tagged when a partner photon of pt> 
500 MeV/c is detected such that the pair falls in an in- 
variant mass window of 108-165 or 500-600 GeV/c^. In 
what follows, all such tagged photon pairs are presumed 
to have been correctly identified as coming from the same 
meson parent. For the rj meson the rate of combinatorial 
pairs is significant and is subtracted using the side-bands. 

Letting the notation tag denote photons which have 
been tagged as coming from decay sources and the no- 
tation miss denote decay photons which could not be 
tagged, this residual decay component can be statisti- 
cally subtracted analogously to Equation [T] according to 
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FIG. 4: The decay probability (normalized by area up to 
Pt = 20 GeV/c) for a tt" to decay into a photon with 5 < 
Pt < 7 GeV/c as a function of with tagging (red) and 
without tagging (black). 
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^inclusive-tag is simply the per-trigger yield of all photons 
remaining in the sample after tagged decay photons are 
removed. The effective for the sample is 



determined to be as large as 5% for the highest p^'^ 
selection. The direct photon contribution to the falsely 
matched sample is assumed to be the same as in the ab- 
sence of photon tagging (i.e., given by Rj) and the total 
size of the direct photon contribution is taken as the sys- 
tematic uncertainty. 



Dmiss -^inclusive -^tag -^7 -^inclusive -^tag 

7 ]\T AT" -1 tag AT ' 

decay -'^tag 1 — ^decay -^^ inclusive 

J7) 

where the tagging efficiency, e^^^l^^, is iVtag/iVdecay Cdecay 
varies from 0.43 in the 5-7 GeV/c bin to 0.53 in the 12-15 
GeV/c bin. 

In order to determine id'ecay' define a decay photon 
probability density p™"'*' in which the decay photon tag- 
ging is performed in the Monte Carlo simulation in the 
same manner as in the data. Figure |4] shows decay prob- 
ability as a function of tt*^ pt for photons in the range 
5 < pI < 7 GeV/c {i.e., P{pl^,pl)dpl)) with and 
without the decay photon tagging applied. The curves 
should be interpreted as the probability (normalized up 
to 7r° Pt = 20 GeV/c) for a tt° to decay into a 5-7 GeV/c 
photon in PHENIX acceptance as a function of 7r° pr. 
The effect of the tagging is most pronounced in this low- 
est decay photon pT bin because the opening angle be- 
tween the decay photons is largest. Photons which pass 
the tagging cut are typically closer to the parent tt" px 
than in the case when no tagging was applied. 

The subtraction procedure defined above assumes that 
all tagged photons were paired correctly. A systematic 
uncertainty is assigned to account for the possibility that 
direct photons are falsely tagged by accidental combi- 
nation with a decay photon. The rate of false tagging is 
estimated from the combinatorial background level deter- 
mined from fits to the invariant mass distributions and 



C. Photon Isolation 

In order to further reduce the decay photon back- 
ground one may impose an isolation requirement, as di- 
rect photons are expected to be largely produced without 
much associated hadronic activity. Such a requirement 
has the additional benefit of partially suppressing the 
fragmentation component. The isolation criterion em- 
ployed in the present analysis is that the sum of the mo- 
menta of charged tracks and the energy of photon clusters 
inside a cone of radius 0.3 centered around the photon 
direction is less 10% of its energy. The cone size is lim- 
ited by the size of the PHENIX aperture which spans 
0.7 units of pseudorapidity. For photons near the edge 
of the detector the isolation cone lies partially outside of 
the PHENIX acceptance. In order to reduce the impact 
of the acceptance on the photon isolation, fiducial cuts 
are applied such that photons are required to be '-^ 0.1 
radians from the edge of the detector in both rj and (p. 

As was the case for decay photon tagging, a residual 
decay component must be statistically subtracted after 
the isolation cuts have been applied. The per-trigger 
yield of isolated direct photon is obtained according to 

■y^iso _ ^ f rpiso y;iso -y^iso \ /o\ 

direct piso Y ^ ' inclusive — tag decay/ \^ ) 

Note that the label iso (niso) denotes photons which 
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are both isolated (non-isolated) and were not removed by 
the tagging cuts. The value of corresponding to the 
photon sample after both types of event-by-event cuts is 



niso 



inclusive 



'^dccay 
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inclusive 



yydccay 
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'^decay 



inclusive 
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decay 
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decay 



(9) 
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The efficiency with which the isolation cut removes 



decay photons, eJJJfj.'i^y, is not known a priori, since an 
unknown fraction of the non-isolated photons are direct. 
In order to evaluate the isolation efficiency we apply the 
isolation cut at the level of the parent mesons and use 
the decay probability functions to map the effect to the 
daughter photon px- For the example of the tt*^, 



1 



Pr^'iPr^Pl) ■ Nr°iPT) 



(10) 



We have implicitly exploited the fact that the tagging 
probability is independent of the isolation requirement. 
e"f°) varies from 0.4 in the 5-7 GeV/c bin to 0.48 in the 
12-15 GeV/c bin. 

The per-trigger yield of isolated (and tagged) tt" decay 
photons can be calculated according to 



viso 



EN^pr'iphPT)Nr-h 

Epmiss /'/ri'^ T\'^ ^ /X/'iso 



(11) 



Azimuthal correlations for inclusive and decay photons 
for 9 < < 12 and 2 <p^"«°'= < 5 GeV/c are shown 
in Fig. [SJ Also shown are the correlations after the tag- 
ging and isolation-f-tagging cuts have been applied. The 
tagging causes a reduction in the near-side yield of in- 
clusive photons confirming that direct production is sub- 
sequently enhanced. The decay photon associated yield 
is also smaller after tagging, due to the fact that the 
remaining photons have a smaller mean px , as can be as- 
certained from the decay probability function shown in 
Fig. m The isolation cut causes a further reduction in 
the yields, driven by the fact that the fraction of the jet 
momentum carried by the trigger is larger for isolated 
photons. Table H] shows the and its effective values 
for the tagged and tagged-|-isolated samples. 



IV. RESULTS 
A. Azimuthal Correlations 

Essential features of dijet or 7-t-jet production are 
evident in azimuthal two particle correlations. Fig- 
ure [H] shows the per-trigger yields of associated charged 
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FIG. 5: Per-trigger yields for inclusive (top) and decay pho- 
tons (bottom) with no event-by-event cuts, with decay photon 
tagging and with decay photon tagging and the isolation cut. 
The systematic uncertainties on the decay photon triggered 
yields are not shown. 



TABLE 1: R-y and its effective values for the photon sample 
with tagging implemented (J?™"^'^) and the sample with tag- 
ging and isolation implemented {Rj°). 



PT [GeV/c] 




nmiss 


r»iso 


5-7 


1.19 ± 0.06 


1.32 ± 0.11 


1.38 ± 0.12 


7-9 


1.33 ± 0.05 


1.67 ± 0.13 


1.92 ± 01.4 


9-12 


1.54 ± 0.05 


2.22 ± 0.18 


2.87 ± 0.22 


12-15 


1.80 ± 0.11 


2.69 ± 0.36 


4.02 ± 0.50 



hadrons as function of A(/) for 7r°, direct photon (using de- 
cay photon tagging) and isolated direct photon triggers. 
The strong near-side correlation for tt'^ triggers is largely 
absent for direct photon triggers as expected from a sam- 
ple dominated by photons produced directly in the hard 
scattering. On the away-side the direct photon triggered 
yields are generally smaller than for tt" triggers. Note 
that the error bars on the direct photon yields, which 
are from R^ and from the decay photon yields, vary bin- 
to-bin depending on the relative values of ^Inclusive and 
^dccay The 8% normalization uncertainty on the charged 
hadron yields is not shown on this plot and all those 
which follow. 
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FIG. 6: Yield per trigger of charged hadrons associated with vr'', direct and isolated direct photons triggers as a function of 
A(/!>. The lines correspond to fits which are described in Section IIV El The constant level, which estimates the contribution 
from the underlying event, is shown as a dashed line (only visible in the lower bins). 



B. Near-side Correlations 

The use of 7 + jet events, whether to calibrate the en- 
ergy of the away-side jet in p -I- p collisions or to study 
its energy loss in nuclear collisions, relies on a detailed 
understanding of photons produced in jet fragmentation, 
which may be considered a background in this regard. 
Evidence of fragmentation photons may be observed indi- 
rectly via near-side correlations of direct photon-hadron 
pairs. Figure [7] shows the ratio of the near-side yield of 
direct photon triggers to that of tt^ triggers for various pt 
selections. The ratio is generally consistent with zero, al- 
beit within fairly larger uncertainties. Such larger uncer- 
tainties on the near-side direct photon per-trigger yields 
are due to the large near-side associated yields for in- 
clusive and decay photon triggers compared to the small 
difference between them. 

The ratio Vdircct/^Tr" is not a direct measurement of the 
relative contribution of fragmentation photons to the di- 
rect photon sample. This quantity is sensitive to both the 
overall contribution from fragmentation photons as well 
as the difference between parton-to-photon and parton- 
to-7r° fragmentation functions and possibly a different 
number of quark vs. gluon jets in the tt^ triggered and 
photon triggered samples. As an illustrative example, if 
one assumes that the photon and tt'' fragmentation func- 
tions are similar in shape, it is reasonable to fit the ratio 



of the per-trigger yields to a constant value as a function 
of 

^assoc rjnj^g rcsults of such fits, shown in Table |lT] show 
that above 7 GeV/c the near-side yield associated with 
direct photons is constrained to be smaller than 15% of 
the TT*^ associated yields. It should be emphasized that 
although such an observation is compatible with a small 
yield of fragmentation photons, it may also reflect a dif- 
ferent shapes of the fragmentation function into photons 
compared to the fragmentation function into neutral pi- 
ous. 



TABLE II: Ratio of constant fits to the near-side yields as a 
function of using isolated direct photon triggers to the 

same quantity using tt" triggers. 



[GeV/c] 


Ratio from fit 


Stat. 


Sys. 


5-7 


0.01 


0.04 


+0.26 - 0.46 


7-9 


-0.03 


0.04 


+0.12 - 0.16 


9-12 


0.04 


0.04 


+0.07 - 0.09 


12-15 


-0.16 


0.07 


+0.14 - 0.17 
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FIG. 7: The ratio of the direct photon to tt*^ near-side {(j) < 
tt/2) associated yields as a function of pj?""'^ . The points are 
placed according to < > for the isolated direct photon 

triggered sample. 



C. Hard Scattering Kinematics using xe and pout 

Due to the effects of hadronization, we do not have 
direct access to the parton kinematics and therefore 
can measure neither the fragmentation functions nor the 
magnitude of the kx effect directly. However, to the ex- 
tent that the LO Compton scattering process (see Fig. [1]) 
dominates, direct photons may be considered to play the 
role of the hard scattered parton. For both the case of iso- 
lated photon and tt'^ triggered correlations we construct a 
simple model to extract the parton-level kinematics from 
the data, as described in the next section. To facilitate 
the interpretation of the data we choose a set of observ- 
ables which are appropriately sensitive to the quantities 
of interest. 

The fragmentation function is not directly measurable 
via two particle correlations since the jet momentum is 
not determined. If the near-side jet pt were fixed by that 
of the trigger particle and the jets were balanced, one 
could measure the fragmentation function of the away- 
side jet by measuring the associated yield as a function 
of the partner px- For the case of dihadron correla- 
tions, however, the of the hard scattering varies with 
both the trigger^and partner px selection, which are both 
jet fragments 16]. Direct photon triggers, on the other 
hand, should provide a nearly mono-energetic sample of 
jets for fixed photon px- The quantity Xh measures the 
transverse momentum imbalance between the trigger and 
associated parton: 



The quantity xe measures the pt balance between the 
trigger and associated particles [SSj] . It is defined in anal- 
ogy to the fragmentation function variable z by substi- 
tuting the the px of the trigger particle for that of the jet 
and taking the projection of the associated particle onto 
the trigger axis in the azimuthal plane: 



Xe 



trig 

_ Pj, 



■Px 



\p 
T 



^assoc I 



trig. 



cos A(f>. 



(13) 



Py'® (which is not the case for jet fragments 



For p^'s 

measured at fixed Py'*), the xe distribution approxi- 
mates the fragmentation function, D{z). Hence, the xe 
distribution for direct photon triggers should scale with 
p^'^ , in the same way as the fragmentation functions 
approximately scale with Q^. 

Beyond LO pQCD the outgoing parton pair may ac- 
quire a net px as depicted in Fig[Sl The pair may acquire 
a net momentum along the both the dijet axis and or- 
thogonal to it. At the hadron level, the kinematics along 
the jet axis are dominated by the effect of jet fragmen- 
tation. The kx effect is therefore best observed by mea- 
suring the momentum of jet fragments in the direction 
orthogonal to the parton pair. Again using the trigger 
particle direction as a proxy for that of the parton, we 
define Pout as a vector transverse to Pj'J"^ of magnitude 



bo 



|p^^''°'=|sinA0. 



(14) 



Pout also contains a contribution from jx which is the 
momentum transverse to the jet axis imparted to the jet 
fragment in the course of the parton showering process. 
Note that the presence of |p|?*'°'^ | in pout implies a depen- 
dence on the fragmentation function of the away-side jet. 
Similarly the longitudinal (along the dijet axis) compo- 
nent of kx can play a role in the xe distributions. Such 
a mixing of longitudinal and transverse effects is an un- 
avoidable consequence of hadronization. 
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FIG. 8: A diagram showing the kinematics underlying the 
measurement of jet correlations between back-to-back parti- 
cles, adapted from p^ . 



xn^pT°'lfx\ (12) 

where the hat symbol denotes partonic level quantities. 
At LO = 1, but may deviate from unity due to the kx 
effect. Given a falling px spectrum, x\i is typically less 
than one due to the trigger bias. 



D. Xe Distributions 

Figure IH] shows the away-side (A0 > 7r/2) charged 
hadron yield per trigger as a function of x e for both tt'^ 
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FIG. 9: Away-side charged hadron yield per tt^' trigger (open 
symbols) and isolated direct photon trigger (closed symbols) 
as a function of xe, which is equivalent to zt in the coUinear 
limit cos{A(l)) = 1. The lines represent fits to the modified 
Hagedorn function described in the text. 



and isolated direct photon triggers. The isolated direct 
photon triggered data are obtained by the same proce- 
dure outlined in Sec. IIIII that was used to obtain the 
A(j) distributions presented in the previous section. The 
isolated direct photon xe distributions are significantly 
steeper than those of tt" . This is to be expected because, 
modulo the smearing, the the x e distribution opposite 
to an isolated direct photon should be the fragmentation 
function of the outgoing parton, predominantly the quark 
from the Compton diagram (Fig. [T]), with a small gluon 
admixture from the annihilation diagram. 

The Xe distributions have been fit using a modified 
Hagedorn function defined by 



dN 

dXE 



N{n - 1) 



1 



Xh (1 



(15) 



where n is the power-law dependence of the inclusive in- 
variant t:^ pt spectrum. This functional form was shown 
to describe the i:^ triggered x e distribution given the fol- 
lowing simplifying assumptions: The hadron is assumed 
to be coUinear with the parton direction, i.e., cos « 1 
in Eq.[T3l the underlying fragmentation functions {D{z)) 
are assumed to take an exponential form and xt is taken 
to be constant as a function of a;^;. The fit was found to 
perform reasonably well for the range below xe < 0.8. 
The results of the fits in this range are shown in Table Hill 
The values arc rather large, particularly at low tt^ pt, 
indicating that the data are not perfectly described by 
the modified Hagedorn function, which is perhaps not 
surprising given the good statistical precision of the tt^ 
triggered data. By further restricting the range of the fit 
to one can obtain values of per degree of freedom much 
closer to unity. For example, a fit to the 5 < p^'® < 7 



GeV/c for 0.3 < Xu < 0.8 yields a per degree of free- 
dom of 17.5/5. However, the extended range of the fit 
allows for a better comparison with the photon triggered 
data, whose poorer statistical precision necessitates fits 
over a larger range. 



TABLE III: Parameters of fits to the xe distributions for tt" 
triggers for xe < 0.8. 



Py"*" [GeV/c] 


N 


Xh 


xVdof 


5-7 


1.25 ± 0.02 


0.71 ± 0.01 


145/7 


7-9 


1.19 ± 0.03 


0.75 ± 0.01 


75/7 


9-12 


1.22 ± 0.06 


0.76 ± 0.03 


18/4 


12-15 


1.33 ± 0.09 


0.75 ± 0.05 


0.05/2 


5-15 


1.24 ± 0.02 


0.72 ± 0.01 


262/26 



For the case of isolated direct photon triggers, Eq. [12] 
will not describe the x e distribution of the away-side jet 
because it is assumed that the trigger particle is also the 
fragment of a jet. However, for a fixed value of n, Eq. [15] 
is a useful quantitative representation of the steepness of 
the X E distribution via the parameter Xh ; a smaller value 
of Xh corresponds to a steeper distribution. Since the pa- 
rameters n and Xh have physical meaning for Tr'^-hadron 
correlations, we keep the same value of n = 8.1 for the 
isolated-7-hadron correlations so as to get a quantitative 
measure of the relative steepness of the isolated-7-hadron 
Xe distribution compared to the Tr^-h distribution. The 
results of the fits, shown in Table lIVl confirm a steeper 
distribution for photon triggers where statistical preci- 
sion allows. 



TABLE IV: Parameters of fits to the xe distributions for 
isolated direct photon triggers. 



[GeV/c] 


N 


Xh 


xVdof 


5-7 


1.35 ± 0.14 


0.60 ± 0.03 


2/7 


7-9 


1.40 ± 0.24 


0.51 ± 0.05 


3/7 


9-12 


0.83 ± 0.18 


0.66 ± 0.10 


0.5/4 


12-15 


0.75 ± 0.24 


0.60 ± 0.15 


1.4/2 


5-15 


1.11 ± 0.07 


0.63 ± 0.02 


44/26 



The independence of the Xh values for 7r° as a function 
of PT suggests that xe scaling should hold for all the 
data combined. This is shown in the left panel of Fig. [lU] 
in which all tt^ pt selections are fit simultaneously. It 
is interesting to note that the failure oi xe scaling in a 
similar measurement (for lower p^, 2-4 GeV/c) at the 
CERN-ISR [3§| , led to the concept of parton transverse 
momentum and Ut- 

For isolated direct photon production, xe scaling is im- 
portant for a more fundamental reason. If the xe distri- 
bution does indeed represent the fragmentation function 
of the opposite parton, then combining all the data (see 
Fig. [TU)) should, apart from NLO effects, give a universal 
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FIG. 10: xe distributions for tt triggers and isolated direct 

photon triggers for all ranges combined, xe is equivalent FIG. 11: ^ = In 1/xe distributions for isolated direct photon 

to zt in the coUinear limit cos(A(?!)) = 1. The dashed line data for all pr ranges combined compared to TASSO mea- 

and solid lines correspond to fits to exponential and modified surements in e+e" collisions at = 14 and 44 GeV. 
Hagedorn (Ea. [T5)l functions, respectively. 



distribution which is a reasonable representation of the 
quark fragmentation function (28j . 

Within the large errors, the xe scaling appears to 
hold. Fits to both Eq. [15] and to a simple exponential 
are shown. The exponential fit (e"''^^) gives the value 
& = 8.2 ± 0.3, with a per degree of freedom of 48/26, 
which is in excellent agreement with the quark fragmen- 
tation function parameterized [l6l . [28l | as a simple expo- 
nential with b — 8.2 for 0.2 < z < 1.0 and inconsistent 
with the value & = 11.4 for the gluon fragmentation func- 
tion. 

Another, recently more popular way [40] to look at 
the fragmentation function is to plot the distribution in 
the MLLA variable [4l| ^ = Inl/z sa \nl/xE which is 
shown in Fig. [TT] The present data compare well to the 
TASSO measurements [i^ in e+ + collisions which 
have been arbitrarily scaled by a factor of 10 to match 
the PHENIX data, which is reasonably consistent with 
the smaller acceptance of the present measurement. This 
again indicates consistency with a quark fragmentation 
function. 



E. Pout Distributions and ^(jpoutP) 

Figure [12] shows the Pout distributions for tt" and iso- 
lated direct photons for the range of 2 < < 10 
GeV/c. The tt" distributions are fit with Gaussian func- 
tions, as well as by Kaplan functions. The Kaplan func- 
tion is of the form C{1 + p^^j/fo)"", where C, n and 
b are free parameters. This function exhibits the same 
limiting behavior at small values of Pout as the Gaussian 
function and transitions to a power-law behavior as pout 
becomes large. The tails of the distributions, above about 
3 GeV/c, clearly exhibit a departure from the Gaussian 



fits. This may signal the transition from a regime domi- 
nated by multiple soft gluon emission to one dominated 
by radiation of a single hard gluon. The isolated direct 
photon data also show an excess above the fit, notably for 
the 7 < p^'^ < 9 GeV/c range. For values of Pout compa- 
rable to p^'^ , the direct photon data are consistent with 
zero yield. This is to be expected on kinematic grounds 
as the momentum of large angle radiation cannot exceed 
the jet momentum, which should be well-approximated 
by the photon momentum. 




FIG. 12: Pout distributions for tt (open symbols) and isolated 
direct photon (filled symbols) triggers. The tt" triggered dis- 
tributions have been fit with both Gaussian (solid lines) and 
Kaplan (dashed lines) functions. Points with bent error bars 
have been shifted slightly for visibility. 

The value of (IfcrP) is determined by measur- 
ing V (boutp)- Rather than obtaining the values of 
\/ (boutp) directly from the pout distributions, fits to the 
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azimuthal correlations are used, following the procedure 
described in [l^, . In this method the contribution of 
the underlying event is parametrized by a A0 indepen- 
dent contribution and the following fit function is used 
to determine the magnitude of Pout from the away-side 
jet width: 



dN, 



real 



1 dN„ 



dA(j) N dA0 



+ C2 ■ 



cWfar 

dAcI) 



(16) 



where 



dA0 



37r/2 
it/2 



cos A(/) 



outy 



(17) 



The underlying event level Co, the near and away-side 
amplitudes Ci and C2 and \/ (|poutp) are free parame- 
ters. The fits to the A</) distributions are shown in Fig. [51 
Using the data for pi^^°^ > 2 GeV/c, we minimize the 
sensitivity to the underlying event whose level is indi- 
cated as a dashed line. Co is determined from the tt'^ 
triggered correlations and treated as fixed for the direct 
photon correlations. Its value was confirmed to be equiv- 
alent for both trigger species in the range 1 < p^^°'^ < 2 
GeV/c where both sets of A(/) distributions could be reli- 
ably fit treating Co as a free parameter. The (|poutp) 
values obtained from the fits are shown in Fig. [131 The 
width of the pout distributions are found to decrease with 
p^'® . At the same value of p^'® , the isolated direct pho- 
ton widths are larger than that of 7r° , which is reasonable 
given that tt'^ triggers on a larger jet momentum. 



V. DISCUSSION AND INTERPRETATION 

A. Leading Order + Gaussian kr Smearing Model 

In order to interpret two-particle correlations of final 
state particles in terms of properties of the hard scattered 
partons, a model using the Born level pQCD cross sec- 
tions and with a Gaussian kr smearing was constructed 
{LO+kr), similar to the model employed in [17]. Using 
this approach, the magnitude of the kr effect is varied 
in the model until the measured values of ^/J}p^mt\^ are 
reproduced. At leading order the differential cross sec- 
tion for back-to-back hadron production from a 2 ^ 2 
scattering process {pi+p2 — > P3+P4) is 



(5, 

V 
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FIG. 13: \/ (Ipout p) for 7r° and isolated direct photon triggers. 



dh 



dxidx2d cos 9* dz'idzi^ 



J2 fa{^i)Mx2 f^jf^ a,Mcose*)D,iz3)Dd{zi). 



a,b,c,d 



(18) 



A Monte- Carlo generator is used to throw flat distri- 
butions of particle pairs in xi,X2,9*, 23 and Z4, where x 
is the fraction of the proton's momentum taken by the 
initial-state parton, 9* is the polar scattering angle, and 
z is the fraction of the final-state parton taken by a frag- 
ment. For pairs which fall in the PHENIX pseudorapidity 
(|?7| < 0.35) interval the right-hand side of Equation [T8l is 
used to weight the contribution of each particle pair. The 
expression consists of the angular component of the LO 
parton scattering cross section, a, for each parton flavor 
combination a+b and the non-perturbative parton distri- 
bution functions (PDF), /, and fragmentation functions 
(FF), D. An equal number of all possible permutations 
of each parton flavor are thrown, neglecting charm and 
heavier quarks. The angular part of the parton scatter- 
ing cross section, i7a,b(cos 9*), contains a numerical factor 
which weights each scattering process appropriately, s is 
the Mandelstam variable representing the square of the 
center-of-mass energy at the partonic level and is related 
to the p + p center-of-mass energy of -^i = 200 GeV by 
s = X1X2S. The amplitudes for dijet processes can be 



found in [4J] . Direct photon- hadron correlations are also 
considered. The 7 + jet amplitudes can be found in [45| . 
In this case, D{z4) is taken to be a J function at Z4 — 1 
and the coupling becomes aga^M- 

The PDFs and FFs are taken from global fit analy- 
ses. The PDFs used in this study are the CTEQ6 set [4I 
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evaluated using NLO evolution which were obtained from 
the Durham HEP database [46|. Several fits to the FFs, 
KKP, AKK and DSS are tested with the Monte Carlo 
generator 113, H^. These parametrizations differ in 
the selection of data sets used in their fits. The KKP set 
relies solely on e+e~ data while the recent AKK and DSS 
fits also employ data from deeply inelastic scattering and 
p + p collisions. The use of these additional data sets en- 
ables the separation of the quark and anti-quark FFs and 
provides much better constraints on the gluon FFs, al- 
though uncertainties for the gluon FFs remain large. De- 
spite the recent progress, the AKK and DSS differ quite 
substantially in a number of observables and discrepan- 
cies exist between fits and certain data sets. Newer data 
sets and observables are required to further constrain the 
FFs. Of particular interest for studies of heavy-ion col- 
lisions is the region of low z where modifications of the 
FF by the dense QCD medium are expected. Whereas 
inclusive hadron cross sections are most sensitive to rel- 
atively large values of z (typically 0.7-0.8), two-particle 
correlations provide access to smaller values of z using 
an asymmetric pr selection of the trigger and partner. 
Direct photons are ideally suited for this purpose due to 
the absence of near-side jet fragmentation effects at LO. 
However, higher order effect such as photon fragmenta- 
tion and soft gluon emission must be constrained. 

The fcr effect is modeled by introducing a Gaussian dis- 
tributed boost of magnitude \p^"\ = V^lkx] randomly 
oriented in the plane transverse to the incoming parton 
pair, i.e. the beam direction. The outgoing parton pair, 
which is back-to-back to leading order, hence acquire an 
acoplanarity and a momentum imbalance. In addition to 
kr the hadrons also acquire some momentum relative to 
their parent parton direction, denoted by jx, from the 
parton showering process. The width of the Jt distribu- 
tion is taken from measurements of Tr'^-hadron correla- 
tions and was determined to be ^ (|jTp) ~ 635 MeV/c 

The Gaussian kx smearing model has the advantage of 
simplicity compared to a more detailed calculation of the 
underlying jet kinematics and subsequent fragmentation. 
Although NLO calculations are available, modeling the 
kx smearing would require care to avoid double-counting 
the contribution at NLO with that of soft gluon radia- 
tion. A self-consistent approach to doing so would require 
a fully resummed calculation for which generator level 
Monte-Carlo simulations are not presently available. Al- 
though a Gaussian kx smearing can be similarly tuned 
in more sophisticated LO models such as PYTHIA, the 
hO+kx model enables us to test various parametriza- 
tions of the fragmentation functions as opposed to using 
a model of the hadronization process. 

It should be emphasized, however, that the Gaussian 
kx smearing model does not take into account the full 
parton kinematics and hence, requires some care to avoid 
unphysical scenarios. The LO cross sections are divergent 
in the forward and/or backward directions and the gluon 
distribution becomes very large at low x. In the absence 



of kx smearing these effects are irrelevant for production 
at midrapidity. However, for a Gaussian distributed kx 
of fixed width there is a finite probability for a parton 
to be scattered at large angle, solely by virtue of receiv- 
ing a large momentum kick from sampling the tail of 
the kx distribution. Due to the largeness of the low x 
gluon distribution and the cross sections at small angle 
these soft partons would dominate the cross section. This 
is clearly an unphysical consequence of the kx smearing 
procedure. The kx boost is intended to simulate gluon 
emission which should clearly be bounded by the momen- 
tum of the parton from which it radiates. This require- 
ment is enforced by imposing the that |A;t| < X\fsj2. The 
calculation was found to be insensitive to the threshold. 



B. Estimating the Magnitude of the kx Effect 

In order to determine the best value of a/ (Ifc^p), 
\J (boutp) is calculated from the model for several dif- 
ferent values of ^ (Ifcrp) and the results are compared 
to data. Figure [14] shows \/ (|poutp) values as well as 
the LO-|-fcT calculation for several values of yj (I^tP) 
and several parametrizations of the FFs. The -k^ trig- 
gered data show that although the LO~|-/ct qualitatively 
reproduces the trend of the data, it does not perfectly re- 
produce the p^^ dependence of \/ (|poutp)- Additional 
effects not included in the LO-t-fc^ smearing model may 
certainly be relevant at this level. At NLO effects such as 
the radiation of a hard gluon may play a role. Alterna- 
tively, soft gluon radiation may not be perfectly described 
by a Gaussian smearing or may depend on the scattering 
process or the momentum exchange of the hard scatter- 
ing. 

The per degree of freedom between data and model 
calculations are show n in Fig . [TSl Direct photon triggers 
give a best value of i\kx\^) ~ 3 GeV/c. The best value 
for TT^ triggers is somewhat larger and depends on the 
choice of FF. This dependence arises from the effect of 
fragmentation on the near-side and the larger fraction of 
gluon jets, for which the parametrizations differ. 

Figure shows \/ (Ifcrp) as function of p^'® as cal- 
culated in the LO-|-fcT model for the best value of the 
input \/ {\kx\'^) as determined from Fig. 1151 The system- 
atic error band indicates the dependence on the choice of 
FF parametrization. The value of \/ (I^Tp) may depend 
on J5^'^ since the trigger requirement may preferentially 
select events based on their kx- For b oth the 7r° and pho- 
ton triggered samples, the value of y/ (IfcrP) depends on 
p^'^ , and reaches values larger than the input, although 
the effect is only statistically significant for the 7r° trig- 
gered data. The values of VxIfcrF) are generally larger 
for TT*' triggers, but within statistical uncertainties on the 
photon triggered sample, the size of the kx effect is of 
comparable magnitude. 
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FIG. 14: \J (IpoutP) compared to LO+fcr calculations for tt' 
triggers (top) and isolated direct photon triggers (bottom) 



C. Sensitivity of xe Distributions to the 
Fragmentation Functions 



Using the value of VvTfcrF) which best matches the 
data, the shape of the xe distribution can be calculated 
in the LO+fcr model, for each set of FFs, and compared 
to the data. Figure [17] (top panel) shows the slope pa- 
rameter, ai/i, of the power- law fits to the xe distributions 
for 7r° triggers shown in Fig. [HI Several calculations are 
shown. The KKP and DSS fits reproduce the shape of the 
data better than the AKK fit, which shows a much harder 
slope. The disparity amongst the calculations seems to 
contradict the claim in [l6| that the x e distribution for tx^ 
triggers is not sensitive to the overall shape of the FF. To 
test this assertion using the LO+fc^ model, an exponen- 
tial function was used for each flavor of FF. The slopes 
of the FFs were varied, using Diz) ex. exp (—8.22:) and 
D{z) (X exp(— 11.4z), to represent the quark and gluon 
FFs, respectively, as was done in fl^. Indeed one finds 
that the calculation is not very sensitive to the change 
in the slope parameter. This exercise does not, however, 
take into account that a change in the shape of the FF 
for an individual parton flavor may change not only the 
admixture of quark and gluon jets, but also the shape of 
the pt distribution of hard scattered partons, which is 
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FIG. 15: between the values of \/ (jpoutp) obtained in data 
and from MC calculation for tt" and isolated direct photon 
triggers. 
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FIG. 16: \/(]fcTpy as a function of Py'^ for tt'' and isolated 
direct photon triggers. Also shown are tt" — ft results using a 
different method [la |. 

taken to be flxed in Since the PDF for the gluon is 
much different than for the quarks, the sensitivity of the 
Xe distribution to the parent parton composition of the 
TfO triggered jets can be tested by removing gluon scat- 
tering processes from the calculation. To illustrate this, 
the slope resulting from using only the q + q ^ q + q 
processes with the DSS FFs is also shown in Fig. [T71 
A signiflcantly harder slope is obtained, verifying that 
the Xe distribution depends on the parton species com- 
position of the triggered jet sample. The effect of the 
kx smearing on the shape of the Xe distribution was 
also investigated. Turning off the smearing results in a 
harder slope for smaller p^'^ , which gradually disappears 
as p^'^ becomes much larger than {kx). 

Figure [T7] (bottom panel) shows the same comparison 
for isolated photon triggers. Here the calculation depends 
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less strongly on the choice of FF parametrization due to 
the smaller contribution of the annihilation subprocess, 
q+q 7 +.9, compared to the predominant g+q — > 7+9. 
In contrast to the 7r° triggered sample, the shape of the 
xe distribution depends rather strongly on the shape of 
the overall FF as demonstrated by varying the slope pa- 
rameter of the exponential parametrization. The magni- 
tude of this effect is reproduced in a gluon jet sample, in 
which only the annihilation process is turned on. The dis- 
tribution for the gluon jet sample is significantly steeper 
than the data, verifying that the sample is dominated by 
quark jet fragmentation. As for the ir^ triggered sam- 
ple, the Monte-Carlo xe distribution becomes harder if 
the kx smearing is turned off. Irrespective of the fit, the 
data in the 5 < p^'* < 7 GeV/c bin are incompatible 
with a model without kx smearing, as a significant yield 
is observed at xe> 1. 
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FIG. 17: Slope parameter Xh of fits to the xe distributions 
shown in Fig. [9] for tt" (top) and direct photon (bottom) trig- 
gers. The slopes are compared to various calculations using 
the LO+kr model as discussed in the text. 



D. Charge Asymmetry 

The dominance of Compton scattering in direct photon 
production implies that the flavor distribution of valence 
quarks in the proton should be reflected in the away- 
side parton. Since there are twice as many up quarks as 
down and the amplitude depends on the electric charge 



of the quark, one expects an asymmetry of 8 : 1 in the 
number of up quark to down quark recoil jets from the 
Compton scattering of a valence quark and a gluon. One 
expects a dilution of this factor for several reasons, for 
example creation of charge pairs in the course of the par- 
ton shower process, corrections from higher order pro- 
cesses such as photon fragmentation and a contribution 
from sea quarks. Nevertheless, a residual charge asym- 
metry should be apparent in the final state hadrons. Fig- 
ure [iH] shows the ratio of positively to negatively charged 
hadrons (i?±) on the away-side of both tt" and isolated 
direct photon triggers as a function of p^^°'^ along with 
calculations using the DSS FFs in the LO+kx model with 
the best value of kr as determined in Section IVB I The 
7r° triggered data show an i?± close to unity. On the 
other hand, an excess of positive charge is evident in the 
isolated direct photon triggered yields. This supports the 
claim that the recoil jet is dominated by quark fragmen- 
tation. 
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FIG. 18: Charge asymmetry of hadrons opposite isolated di- 
rect photon and 7r° triggers as a function of a:_B and the cor- 
responding calculations in the LO+kr model using the DSS 
FFs. 



VI. CONCLUSIONS 

A detailed understanding of jet fragmentation in p + p 
collisions is a prerequisite for studies of possible medium 
modifications to the fragmentation functions by the QCD 
medium at RHIC. To this end, this study provides base- 
line measurements of two particle correlations from which 
to quantify these effects. Direct photon triggered cor- 
relations were obtained by a combination of event-by- 
event identification, in the form of decay photon tagging 
and isolation cuts, and a subtraction of the residual de- 
cay photon associated component. The measurement of 
near-side hadron production associated with inclusive di- 
rect photon (i.e., without isolation cuts) triggers shows 
no evidence for a large contribution from dijet processes. 
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as might be expected from a large fragmentation pho- 
ton component. The yield opposite isolated direct pho- 
tons was found to be smaller than for tt^, consistent with 
the expectation of a smaller away-side jet momentum for 
fixed py'^ . Furthermore, the isolated direct photon asso- 
ciated yields were steeper than for tt" and demonstrated 
Xe scaling as one would expect if the distributions closely 
resemble the underlying charged hadron fragmentation 
function. As a function of Pout, the away-side yields for 
7r° triggers show a Gaussian-like behavior at small val- 
ues of Pout, whereas a harder, power-law like component 
emerges at large pout • The tail component is interpreted 
to be due the emission of a hard gluon. The isolated di- 
rect photon triggered distributions also appear to show 
evidence for a hard tail; however the data are consistent 
with zero yield when Pout > Pj''^ > corresponding to the 
kinematic limit for LO photon production. 

The results were further interpreted at the parton 
level using a simple model of LO pQCD incorporating 
a phonomcnologically-motivated Gaussian kx smearing. 
The hadron yields opposite isolated direct photons are 
shown to be directly sensitive to the fragmentation func- 
tion of the away-side parton. In contrast, hadron trig- 
gered jets are sensitive only indirectly, due to the contri- 
bution of multiple sub-processes in the initial state, the 
relative contribution of which is sensitive to the gluon 
FFs. Furthermore, the shape of the distributions are 
shown to be compatible with the Compton scattering 
process q + g ^ 9 + 7- The dominance of the Compton 
scattering process is further reinforced by the positive 
charge asymmetry observed opposite isolated direct pho- 
ton triggers. Finally, the d irect ph oton data are shown 
to be compatible with a -\/ (IfcrP) of similar magnitude 
to that required by the tt'^ triggered data, within un- 
certainties. Such a large momentum imbalance between 
the photon and the recoil jet is significant for studies of 
photon tagged jets in nuclear collisions in the kinematic 
regime currently accessible at RHIC. 
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